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ON JUMP-DIFFUSION PROCESSES WITH REGIME 
SWITCHING: MARTINGALE APPROACH 

ANTONIO DI CRESCENZO AND NIKITA RATANOV 


Abstract. We study jump-diffusion processes with parameters switching at 
random times. Being motivated by possible applications, we characterise 
equivalent martingale measures for these processes by means of the relative 
entropy. The minimal entropy approach is also developed. It is shown that in 
contrast to the case of Levy processes, for this model an Esscher transformation 
does not produce the minimal relative entropy. 


1. Introduction 

We investigate some basic properties of the jump-diffusion processes 
X(t) = T c (t) + N h (t) + W a (t), t> 0, 

with time-dependent deterministic driving parameters switching simultaneously at 
random times. Here W a denotes the Wiener part, defined by the stochastic integral 
(w.r.t. Brownian motion B) of the process which is formed by switching at random 
times of the deterministic diffusion coefficients cr* = c7j(t), i £ D := {1,..., d}, d > 
2; by N h is denoted the jump part, i.e. the stochastic integral (w.r.t. process 
N = N(t ) counting number of the regime switchings) applied to the switching 
functions hi = hi(t ), i € D, and T c is path-by-path integral in time t of switching 
velocity regimes c* = Ci(t), i £ D (see the detailed definitions in Section 2). In 
the case of d = 2 and exponentially distributed inter-switching time intervals such 
processes are called telegraph-jump-diffusion processes, Ratanov (2010), or Markov 
modulated jump-diffusion. 

In this paper the random inter-switching time intervals are assumed to be in¬ 
dependent and arbitrary distributed. In general, such a process is not Markovian, 
and it is not a Levy process as well. We study these processes from the martingale 
point of view, including Girsanov’s measure transform. 

Similar models without a diffusion component were considered first in Melnikov 
and Ratanov (2007) and more recently have been analysed in detail by Di Crescenzo 
and Martinucci (2013); Di Crescenzo et al (2013), Ratanov (2015, 2013, 2014a), 
see also the particular cases in Di Crescenzo (2001); Di Crescenzo and Martinucci 
(2010). The model with missing jump component is presented in Di Crescenzo 
ct al (2014); Di Crescenzo and Zacks (2015). The processes with random driving 
parameters are studied in Ratanov (2013). The recent paper Ratanov (2014b) is 
related to the model of random switching intensities. 
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This setting is widely used for applications, see e.g. Weiss (1994). The mar¬ 
tingale approach developed in this paper is motivated by financial modelling, see 
Runggaldier (2004) for jump-diffusion model. See also Ratanov (2007) for jump- 
telegraph model (and a more detailed presentation in Kolesnik and Ratanov (2013)). 

The Markov modulated jump-diffusion model of asset pricing (with additive 
jumps superimposed on the diffusion) has been studied before, see Guo (2001); 
Ratanov (2010). This model for a single risky asset possesses infinitely many mar¬ 
tingale measures, and thus the market is incomplete. The model can be completed 
by adding a further asset; for the jump-diffusion model see Runggaldier (2004), and 
for the telegraph-jump-diffusion model see Ratanov (2010). 

In this paper we explore another approach. We describe the set of equivalent mar¬ 
tingale measures and determine the Follmer-Schweizer minimal probability mea¬ 
sure (so called the minimal entropy martingale measure (MEMM)), Follmer and 
Schwcizer (1990). In his seminal paper Fritelli (2000) Fritelli has showed the equiv¬ 
alence between maximisation of expected exponential utility and the minimisation 
of the relative entropy. Then, by this approach the models based on Levy processes 
have been studied in Fujiwara and Miyahara (2003). 

Observe that for Levy processes and for regime switching diffusions the usual 
technique is based on an Esscher transform, which produces the MEMM, see Fu¬ 
jiwara and Miyahara (2003); Esche and Schweizer (2005) and Elliott et al (2005, 
2007). 

In our model this method does not work. The Esscher transform under regime 
switching does not affect the switching intensities. In the Levy model hijiwara and 
Miyahara (2003) this contradicts the minimal entropy condition (if the process is 
not already a martingale), see (5.4). In the case of the regime switching model 
Elliott et al (2005) some additional entropy given by the jumps embedded into this 
model can be reduced by more flexible measure transformation, see Section 5. 

Moreover, in contrast with Levy model and with the regime switching diffusions 
without jumps our model has the following important feature: the entropy minimum 
as well as calibration of MEMM depend on the time horizon under consideration. 

The paper is organised as follows. Section 2 contains the definition and the main 
properties of the regime switching jump-diffusion processes with arbitrary distri¬ 
butions of inter-switching time intervals. We construct Girsanov’s transformation 
in Section 3. Then, we define the entropy and derive the corresponding Volterra 
equations. In Section 4 we describe the set of equivalent martingale measures for 
the regime switching jump-diffusion processes. The minimal entropy equivalent 
martingale measures are studied in Section 5 for the case of constant parameters. 

Some applications including financial modelling will be presented elsewhere later. 


2. Generalised telegraph-jump-diffusion processes. Distributions and 

EXPECTATIONS 

Let (fl,.F, P) be a complete probability space with the given right-continuous 
filtration F t , t > 0, satisfying the usual hypotheses, Jeanblanc and Chesney (2009). 
We start with a d-state J^-adapted semi-Markov (see Jacobsen (2006)) random 
process e, e(t) £ D, t > 0. The switchings occur at random times r n , n > 0, r 0 = 
0, and process e is right-continuous with left-hand limits. 

Let N = N(t) be the counting process, N(t) = max{n | T n < t}, t > 0. 
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2.1. The definition of telegraph-jump-diffusion process. For the set of de¬ 
terministic measurable functions Zi = Zi(t ), t > 0, i £ D, we construct first the 
piecewise deterministic random process combining Zi, i £ D, by means of the 
switching process e: 


z \t) = J2 Z e(T n -i)(t-T n -i)l An (t), t> 0. (2.1) 

n= 1 

Here A„ := [r„_i, r n ), n > 1, and Ia(-) is the indicator function. Process z' 
starts from the origin at the switching time tq = 0; at each further switching time 
r n , n > 1, the process is renewed. 

Second, consider the integrals of z t of the following three types: 

(1) the generalised d-state telegraph process (path-by-path integral) 

T z (t) = f z\u)du, t> 0; (2.2) 

J o 

(2) the pure jump process (integral w.r.t. counting process N) 

N z [t) = [ z\u)dN(u), t> 0, (2.3) 

Jo 

(3) the Wiener process (Ito integral) 

W z {t) = [ z\u)dB(u), t>0, (2.4) 

Jo 

where B = B(t), t > 0, is an ^-adapted Brownian motion, independent of 
e. Note that W z = W z (t), t > 0, is a Gaussian (J 7 t ,P)-martingale. 

Let Ci = Ci(t), hi = hi(t) and cr = CTj(t), t > 0, i £ D , be deterministic 
measurable functions. We assume that functions are locally integrable, and cr 
are locally square integrable, 

(Ti(u) 2 du < oo, t > 0, i £ D. 

In this paper we analyse the jump-diffusion process with switching regimes X = 
X(t), t > 0, of the following form 

X(t) :=T c {t) + N h (t) + W a (t), t> 0, (2.5) 

with the components T c , N h and W a which are defined by (2.2), (2.3) and (2.4) 
respectively. Process X satisfies the stochastic equation 



dX(f) = c EiTNm) {t-T N{t) )dt+h eiTNm) (T N{t) )dN(t)+a £{TNm) (t-T N{t) )dB(t), t > 0. 
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Equivalently, processes T c , N h , W a can be expressed by summing up the paths 
between the consequent switching instants r n : 

,t mt) 

T c {t) = / c\u)du = Y + l s (r NW )(t ~ T N(t)), (2-6) 

n= 1 

,t JVC*) 

N h (t) = / h\u)dN(u) = Y h e(rn-^Tn), (2.7) 

n—1 

,t JVC*) 

w a {t) = / a\u)dB{u) = Y ^(Tn-O^n) + w e (T NW )(t - Tjv(t)), ( 2 . 8 ) 

n—1 

where the following notations are used 

li(t) = f Ci(u)du, Wi(t) = f ai(u)dB(u), t > 0, i G D. (2.9) 
Jo Jo 

Note that h(t), i £ D, are deterministic and the variables Wi(t), i £ D, are zero- 
mean normally distributed with the variances Ei(t ) 2 := f Q cri(u) 2 du, t > 0, i £ D. 

We say that the regime switching jump-diffusion process A' defined by (2.5)- 
(2.9) is characterised by the triplet (c, h, a) with distributions of inter-switching 
time intervals T n which are determined by the hazard rate functions 7 ,y = 7 ij(t), 
i,j £ D , (see the definition in ( 2 . 11 )). 

We apply the notations Xi and T[, . Wf, if the initial state i £ D of the 

underlying process e is given, e( 0 ) = i. 

Further, we will need the following explicit expression for the stochastic expo¬ 
nential £AX) of X =T C + N h + W a . It is known, see e.g. Jeanblanc and Chesney 
(2009), that 

N t 

£t(X) = exp (r c -° 2 / 2 {t) + W°(t)) JJ(1 + h. e(Tn) {T n )) = exp (Y(t)), (2.10) 

n =1 

where Y(t) = T c ~ a J 2 (t) + N ln( - 1 + h '>(t)+W a (t), t > 0. Here T c ~ a / 2 is the telegraph 
process (2.6) with the velocities c, — of/2 instead of Cj, i £ D, and N ln A+ h ) i s the 
pure jump process (2.7) with the jump values ln(l + hi), instead of hi, i £ D, 
switching at random times T n , n > 1 . 

2.2. Semi-Markov process e = e(t). In order to state the distribution of A1(f) 
we introduce conditions on the driving processes e and {r n }n>i. 

Denote by Fij = Fij(t), t > 0, i,j£ D, the transition probabilities of the form 

Jacobsen (2006, (3.17)): 

Fij(t) = P{ti > t, e(n) = j I e(0) = *}, t > 0, i,j £ D. 

Let Fij(t) > 0, t > 0, i,j £ D. Consider the hazard functions Ey, 

r b(*) = _ln - F b(*)> t>0,i,j£D. 

— d ~Fi ■ 

Let functions be differentiable, u ( t) = —fijit), t > 0, i,j £ D. Thus the 
hazard functions 1 /j = Tij(t) are expressed by 
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Here 


7 ij(u) ■■= = 


fij fa) 


Fij (u) ’ 


u > 0, i,j £ D, 


( 2 . 11 ) 


are the hazard rate functions, and fij , i, j £ D 1 are the density functions of the 
interarrival times. We assume the non-exploding condition to be hold: 


pOO 

/ 7 jj(u)du =+ 00 , i,j 
Jo 


£ D. 


( 2 . 12 ) 


Note that 

Fij(t) = exp ( - / ^fij(u)du ) , fij(t) = jij (t) exp ( - / 7 ^(u)d u ) , t>0, i,j G D. 


'0 


/ o 


(2.13) 


The survival function of the first switching time n is given by 

rt 


Fi(t) := P{n > £1 e( 0 ) = *} = P (£) = exp J 7 j(u)du) , t > 0 , i £ D, 

(2.14) 


3£D\{i} 


where 


Furthermore, let 


7* = 53 7u> 

ien\{z} 


fi(t) = -^T(t) = 7 i(t) exp 7 i(u)d, £ > 0 , 


* G D. 


(2.15) 


(2.16) 


Due to (2.13) iV is an inhomogeneous Poisson process with switchings at r n , n > 
1, and with the instantaneous intensities 7 ,j(f), £ > 0, i,j £ D. 

If the process e is observed beginning from the time s, tq < s < n, the cor¬ 
responding conditional distributions can be described by the following conditional 
survival functions, 


Fij(t | s) =P(n > £, e(n) = j | n > s, e(0) = i) = -= 


Fij(t) 


Fi(s) ’ 


(2.17) 


with the density functions 




f i j(t\s) = -—F ij (t\s) = = L —, 0 < s < £, i,j£D. 

Moreover, let 

Fj(£ | s) := P(ri > £ j n > s, e(0) = i) =exp J 7 j(u)du 

fi(t | s) = 7*(£) exp 7i(«)du 

see (2.14). 

Notice that 


i £ D, 
(2.18) 


Fij (£ | 0) =Fij(t), Fi(t | 0) = Fi(£), 

fij(t | 0) =fij(t), fi(t | 0) = /*(£), 


£ > 0, £ D. 
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Assume the inter-switching time intervals T n = A r n = r n — r„_i, n > 1, to be 
independent, and process e to be renewal in the following sense: n > 1, 

P{T n > t | T n > s, e(r n _i) = ij = Fi(t | s), t > s > 0, iSd, 

see (2.17)-(2.18). 


2.3. The distribution and expectation of X(t). Our further analysis is based 
on the following observation. 

Let t i = t be the first switching time, where * € D is the fixed initial state of 
the process e, e(0) = i. Owing to the renewal character of the counting process N, 
see Cox (1962), we have the following equalities in distribution: 


Xi(t) — ( li(t ) + Wj(t))l Tl >t + 


k{n) + Wi(n) + hi(r\) + X e .( Tl )(t 





t > 0, ie D, 


(2.i9) 

where X is the regime switching jump-diffusion process independent of X which 
starts at time ti. Here 1a is the indicator of event A. 

Let 


Pi(x,t) := P(Xj(t) £ da:)/dx, t>0, 
Pi(x,t | s) := P(Xj(t) £ dx | Ni(s) = 0)/dx, t> s, 


x £ (—oo, oo), i £ D, 


be the density functions of Xi(t), i £ D. Note that Pi(x,t \ 0) =pi(x,t), i £ D. 
In these terms equalities (2.19) take the following form 


Pi(x,t | s) =ipi(x-li(t),t)Fi(t | s) 

+ ^2 / (/ Pj(x-k(u)-hi{u)-y,t-u)ipi{y,u)dy) fij(u \ s)dh 

t > s > 0, x £ (—oo, oo), i £ D. 

( 2 . 20 ) 

Here ipi = ipi(x, t) is the density function of the Gaussian random variable Wi(t ), t > 

0, i £ D, 


where Si(t) 2 := / n * ai(u) 2 du, t > 0, i £ D. 

In the Markov case of two-state processes, d = 2, with constant parameters 
Ci,hi,cri, i £ {1,2}, and with exponentially distributed inter-switching times T„, 
the distributions of X\ (t) and X 2 (t) have been analysed in detail in Ratanov (2010). 
Let us study the expectations 


■= E[Xj(f)], t > 0, i £ D, 

and 

Hi{t | s) := E[Xj(f) | Ni(s) = 0], t > 0, i £ D, 
of the d-state process X = X(t), t > 0. 

Note that E[TT CT (t)] = 0, Vf > 0. Moreover, for t £ [0, s], pi(t \ s ) = k(t), i £ D. 







ON JUMP-DIFFUSION PROCESSES WITH REGIME SWITCHING 


7 


To characterise Pi(t | s), i £ D , we use the following notations: 

ai(t) := f [a{u)Fi{u) + hi(u)fi(u)\ dit, t > 0, 

Jo 


i £ D. 


di(t | s) := Zi(a) + / [ci( , u)T 1 j(n | s) + hi(u)fi(u | s)]dit, t > s > 0 , 


Here fi(u), F *(u) and fi(u \ s), Fi(u | s), i G D, are defined in (2.14), (2.16) and 
(2.18). 

Proposition 2.1. The expectations Hi(t), t > 0, i G D, satisfy the following 
Volterra system of integral equations , 

M»(*) = a»(*)+ y: 

If functions Pi(t), t > 0, i £ D, solve system (2.21), then the conditional expec¬ 
tations Hi(t | s) are given by 

Hi(t | s) = ai(t | s) + ^ 


Pj(t — u)fij(u | s)du, t > s, i G D. (2.22) 


/ Hj(t- u)fij(u)du, 
Jo 


t> 0, i £ D. 


( 2 . 21 ) 


Proof. By applying (2.20) (with s = 0) to 

/ OO 

zpi(x,f)da; 

-OO 

one can get 


rt / /»oo /*oo 

Pi{t) =li{t)Fi(t) + / (/ ‘ l l , i(y> u )dy xpj(x - k(u) - hi(u) - y,t - u)dx) fij 

=li(t)Fi(t)+ V] [ [jj,j(t-u) + li(u) + hi(u)]fij(u)du. 

™ r. 1 JO 


Integrating by parts we have 


Pi(t) = k(t)Fi(t)+ V] f (fij(t - u) + hi{u)) fij(u)du- lift)F^(t)+ f Ci(u)Fij(u)di 

™ r , OO JO 




which gives (2.21). The proof of (2.22) is similar. □ 

In the case d = 2 equations (2.21) are derived e.g. in equations (3.2)-(3.3) of 

Ratanov (2013). 

Corollary 2.2. The identities 

Hi(t) =0, t > 0, and pi{t | s) = U(s ), t>s, i £ D, 
hold if and only if 


7 i{t)hi(t) +a(t) = 0, t > 0, i £ D, (2.23) 

where 7 $ = 7 j(t), i > 0, i £ D, are the hazard rate functions , which are defined by 
(2.11), (2.15). 


(u)du 
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Proof. Notice that systems (2.21) and (2.22) have the trivial solutions Pi(t) = 0, 
t > 0, i € D and p.i(t | s) = Ufs), t > s, i £ D, respectively, if and only if cii(t) = 0 
and cn{t | s) = lo(s), i € D. These equalities hold, when 

Ci{u)Fi(u) + hi (u)fi(u) = 0, 

_ «>s>0,ieD, 

Ci(u)Fi{u | s ) + hi(u)fi(u | s) = 0, 

which is equivalent to (2.23), due to (2.11)-(2.13) and (2.17). □ 

Remark 2.3. Condition (2.23) has the sense of Doob-Meyer decomposition. These 
type of conditions for the jump-telegraph processes appears first in Ratanov (2007, 
Theorem 1) in the case of constant deterministic parameters c,h,j ( see also 
Kolesnik and Ratanov (2013)). In this case condition (2.23) is intuitively evident. 
It means that the displacement performed by the telegraph process during a time- 
period r equal to the mean-switching-time is identical to the jump’s size performed 
in the opposite direction. 

This intuitively explains why this is a martingale condition. 

3. Girsanov’s transformation 

In this section we analyse the problems which are important for applications, e.g. 
for financial modelling. First, we describe all possible martingales in our setting. 
Then, we derive a generalisation of Girsanov’s Theorem. 

3.1. Martingale’s characterisation. Since the diffusion part W a = W a (t) = 
J* (u)dB(u) is already a P-martingale, it is sufficient to investigate the process 
X = T c {t) + N h (t), t > 0. 

Theorem 3.1. Let X — X(t), t > 0, be the process with the parameters (a, hi ), i £ 
D 1 switching at random times r„, n > 0. Let the inter-switching times T n = r n — 
t„_i, n > 1, be distributed with hazard rate functions 7 i = 7i(t), t> 0, i € H, see 
(2.11), (2.15). 

Process X is a martingale if and only if the equalities in (2.23) are fulfilled. 

Proof. If A' is a martingale, then no(t) = (t) = 0, which is equivalent to (2.23), 

see Proposition 2.1. 

Conversely, it is known, see Jeanblanc and Rutkowski (2002), Proposition 2.13, 
that the compensated jump process N h — T lh is a martingale. Therefore, if iden¬ 
tities (2.23) hold, then —T c is the compensator of N h and the sum T c + N h is a 
martingale. □ 

Notice that if jumps vanish, hi = 0, i € D, process A is a martingale only in 
the trivial case: a = 0, i € D, and thus A = 0, see (2.23). 

Corollary 3.2 (Ratanov (2010)). Let X = X(t), t > 0, be the jump-diffusion 
process with switching constant parameters Ci , hi and <Ti, i £ D. Let hi ^ 0, i £ D. 

Process X is a martingale if and only if Ci/hi <0, * £ D, and the distributions 
of the inter-switching times T n are exponential, Exp(A,;), with parameters Aj, A, = 
—Ci/hi >0, i £ D. In this case the underlying £ is a Markov process. 

Remark 3.3. In the paper by Di Crescenzo et al Di Crescenzo et al (2014) the 
generalised 2-state geometric telegraph-diffusion process S = S(t) with constant 
parameters Co, ci and a is studied, 

S(t) = s 0 exp [T(t) + crB(t )], 
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where B is a standard Brownian motion and the inter-switching times are indepen¬ 
dent and arbitrarily distributed. The authors expected that the process S = Sit) 
can be transformed in a martingale by superimposing of a jump component. This 
expectation is not justified. 

The process S(t)/so is the stochastic exponential of X = T(t ) + crB{t) + a 2 t,/2. 
After the inclusion of a jump component with constant jump amplitudes hi,h ,2 > 


— 1, such that 


Cj + cr 2 /2 
hi 


< 0, i £ D = {1, 2}, processes X and S become martingales 


only in the standard case of exponentially distributed inter-arrival times, Exp(Aj), 
with constant intensities 


K 


Ci + < t 2 /2 
hi 


i€ D = {1,2}, 


(see Corollary 3.2). 


3.2. Girsanov’s Theorem. The problem of existence and uniqueness of an equiv¬ 
alent martingale measure is extremely significant for applications, especially in the 
theory of financial derivatives. It is important to understand how the equivalent 
martingale measures can be constructed if such a measure exists. 

Let e(t) G D, t > 0, be the switching process on the filtered probability space 
(f2, X, JFt,W) governed by the hazard rate functions 7 i = 7f(t) = fi(t)/Fi(t ), t > 
0,i G D, of the inter-switclring times, see (2.15), (2.11), satisfying the non-exploding 
condition (2.12). 

Let c* and h*, i G D, be measurable functions satisfying the martingale condition 

(2.23), 

7? (t)h*(t) + c*(t) = 0 t > 0, i G D. (3.1) 

We assume c*, i G D, to be locally integrable and h*(t) > —1, t > 0, i G D. Thus, 

c*(t) < 7?(t), Vt > 0, i G D. (3.2) 

Furthermore, let 

(c*(u) — Ji{u)) du = — 00, i G D. (3.3) 

Consider the jump-diffusion martingale A* = T*(t) + N* (t) + W* (t) with regime 
switching, where 

N(t) 

= E + l* ei r NW )<t - TN(t)), 

n =1 
N(t) 

7V*(t)=^/ l *( Tn _ i) (T„), (3.4) 

n =1 
N(t) 

W*(t ) = E <(.r n -,){Tn) + W* e{TNw) {t - T N{t) ). 

n =1 

Here, see (2.9), 

c*{u)du, w*(t)=[ a*(u)dB(u), 

Jo 

where <j*, i E D, are locally square integrable functions. 
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Let Z = Z(t) = £t{X*) be the stochastic exponential of X*. By (2.10) 

N t 

Z(t ) = £t(X*) = exp (T c *- CT * 2 / 2 (t) + W°\t)) na + K(r n ^)(Tn)) 


71=1 


(3.5) 


:exp ( t c *- CT * 2 / 2 (t) + iv ln( 1+/l *) + TT CT *(t)) . 


Here T c ~ <T 2//2 is the d-state generalised telegraph process (2.6) with the velocity 
regimes c* — (<t*) 2 /2 instead of c*, and 7V ln ^ 1+,^ *^ is the pure jump process (2.7) 
with the jump values ln(l + h*) instead of hi, i £ D. 

Theorem 3.4 (Girsanov’s Theorem). Assume that conditions (3.1)-(3.3) hold. 

Let measure Q be equivalent to P under the fixed time horizon t, t > 0, with the 
density 

= (3-6) 

Under the measure Q : 

(a) the inter-arrival times {T n = r n — T ra _i} ra >i are independent and dis¬ 
tributed with the survival functions 

t>0,i£D. (3.7) 

The hazard rate functions 7 ® of these distributions are given by 

7 ?{t) = 7r(i) - c* ( t ) = (1 + h* (t)) 7 f (: t ), t > 0, * G D, (3.8) 

and the non-exploding condition 


pOO 

/ 7?( M )dM = + 00 , i £ D, 

Jo 


(3.9) 


holds; 

(b) the process B(t ) = B{t) — L*{t ) is the standard Q -Brownian motion, 
where L*(t), t > 0, is the generalised telegraph process with switching ve¬ 
locities a*, i £ D, i.e. 

N (t) r T n 


L*(t) '■= T a * (t) = E / "<(r„_ 1 )(«)dw+ / 
J 0 Jo 


t — T N (t) 


£ ( r lV(t)) 


(u)du; 


(c) the ¥-jump-diffusion process X with switching regimes, which is defined 
by (2.5), 

X = T c [t) + N h {t) + W a (t), t > 0, 
under measure Q is still a jump-diffusion process 

X = T c+<Tcr * (i) + N h (t) + W a (t), t > 0, 

characterised by (c + aa*,h,a), with the hazard rate functions 7 ® of the 
inter-switching times {T n } n > 1 , determined by (3.8). Here W a is the sto¬ 
chastic integral (2.8) based on the Q-Brownian motion B. 

Proof. By definition, see (3.5)-(3.6), we have 

7f(t) =Q{n > t | e(0) = *} = E P {Z(t)l {T1>t} | e(0) = t} 


= exp 


c*{u) - {a*)\u)/2]du E [exp (w*(t))) P(ti > t \ e(0) = i). 
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Owing to E [exp(ui*(t))] = exp ^4 f*(cr*) 2 (u)duj we obtain (3.7). 
Note that by (3.7) and (3.3) 


exp 



= F ? (t) = exp (l*(t))Ft(t) = exp 



^(u))duj , 
(3.10) 


t > 0, * € D. 


Hence, 7 ® = 7 ? — c*. Since by (3.1) c* = —h* yf, and (3.8) is completely proved. 
The non-exploding condition (3.9) follows from (3.3). 

The part (b) of the theorem follows from the classical Girsanov’s Theorem, see 
e.g. Jeanblanc and Chesney (2009), Proposition 1.7.3.1. 

The part (c) holds by the following observation. The Wiener part W° of process 
X under measure Q becomes W a (t) = W a (t) + T CTCT *(f), see part (b). Here W a is 
the Q-Wiener process, i.e. the Ito integral w.r.t. Q-Brownian motion B, and T aa 
is the Q-telegraph process which is driven by the subsequently switching velocities 
c Ti(t)cr*(t ), i G D. 

Therefore, under measure Q the process X is still the jump-diffusion process 
with switching regimes, 


X(t) = T c+GrT * {t) + N h (t) + W"(t), t > 0, 


characterised by (c + a a *, h, a). The theorem is proved. 


□ 


3.3. Relative entropy. Let P and Q be two equivalent measures. Under the time 
horizon t, t > 0, the relative entropy of Q w.r.t. P is defined by the set of functions 
Hi(t), t > 0, i £ D : 


Hi(t) := E(£ 


dP^ ; 


e( 0 ) = i 


= Ep 


dP (i) ” dP (t) 


e( 0 ) = i 


, (3.11) 


see Frit ;lli (2000). Here the Radon-Nikodym derivative ~^{t) = £ t (X*) is pre- 

dP 

sented by (3.5)-(3.6). 

Theorem 3.5. Let conditions (3.1)-(3.3) hold. 

The relative entropy functions Hi are expressed by 


Hi(t) = Eq [if +(<T * ) 2 / 2 (t) + Nl n(1+h '\t) 


and satisfy the system of the integral equations 


t> 0, i € D, 


(3.12) 


Hi(t) = ai(t) + ^ f Hj(t — u)fij(u)du, t > 0, i € D. (3.13) 

,,i0 


j€D\{i} 

where functions ai are defined by 


a i{t) = / bi(u)F,i (u)du, t > 0, i £ D. 

Jo 


(3.14) 
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Here 


F?(u) =exp 7 ^(u')du^ , 

b i(u) = 7 * P (w) - 7 ®(u) + + \a*(u) 2 , 


4>i{u) = 


7 “ (it) In 


77C») 

7f(«) 


if 7 i O) ± 0, 

0, i/7z P (w) = 0, 


(3.15) 


it > 0, i £ D. 


Proof. Owing to (3.5) 




ln dF (i) 


e(0) = i 


= Ein 


T i c * _ ( CT *) 2 / 2 (t) + Nl n(1+h *\t) + Wf*(t) 


where the alternating tendencies T i 


c *-( CT *) 2 /2 


(3.16) 

, the jump process N^ n ( 1+h *^ and the 
Wiener process Wf* are defined by (3.4) (with c* — (cr*) 2 /2 instead of c* and 
ln(l + h*) instead of h*, i £ D). 


By Theorem 3.4, part (c), the process T- 


,c*-(<t*)V2 


(i) + N^ 1+h '\t) + Wf*{t) 


nC* + (cr*) /2^ jyln(l+/i* 

the stochastic integral w.r.t. the Q-Brownian motion B. Therefore, E(r 
0, t > 0, and 

Hi(t) =E<(; 

=Er 


under measure Q becomes Tf ’ 'it) + N. n '' l ~ tn ‘ 1 ( t ) + Wf' it), where Wf’ it) is 

r wf(t ) 


t(~ +{ a * )2/2 it) + Nl n(1+h *\t ) + Wf{t) 


t: 


■ c *+(< t *) 2 /2 


(t)+JV] n(1+ ' l * ) (t) 


which gives (3.12). 

Equations (3.13)-(3.15) follow from Proposition 2.1 and Theorem 3.4, see (3.8). 

□ 

It is easy to see that functions hi defined by (3.15) are non-negative, 6 j(u) > 
0, u > 0, i £ D. Hence, functions a* defined by (3.14) are also non-negative, 
a,;(t) >0, t > 0, i £ D. 

Remark 3.6. By applying the Laplace transform / —>• /(s) = J 0 °° e~ st fit)dt to 
(3.13) one can obtain the system: 

Hiis)=aiis)+ ^2 fijis)Hjis), s > 0, i £ D. 

j£D\{i} 

if the transformations cq(s), i £ D, exist. The above system yields the unique 
solution. 

For example, if d, = 2, and b 1 , &2 (see (3.15)) are constants; if the alternating 
distributions of inter-arrival times are exponential, 7 ®(f) = A* = const, i £ {1,2}, 
therefore in this simple case 

B Hi 
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where 


-^l — 


AK&1-&2) 

M+A£) 2 ’ 


Ao = 


a ;(b 2 -M 

(A^ + A2) 2 : 


B = 


AJ + A 2 


(3.17) 


This corresponds to the following explicit solution of (3.13): the relative entropy 
functions Hi(t), H 2 (t), t > 0, are expressed by 


H 1 {t) = H 1 {t- At, AJ) =Bt + Ai 
H 2 (t) = H 2 (t- At, A|) =Bt + A 2 


J _ g — (AJ+AjJjt 
J_ _ g — (AJ + AJ)t 


(3.18) 


4. Equivalent martingale measure 

Consider the jump-diffusion process X = T c (f) + H h (t ) + H /fT (f), t > 0, with 
the switching hazard rate functions yf, i £ D 7 of inter-arrival times T n , n > 1, see 
the definitions in ( 2 . 6 )-( 2 . 8 ). 

Let the equivalent measure Q be defined by the Radon-Nikodym density Z(t ) = 
t > 0, see (3.5)-(3.6). Let driving parameters c*, h*, h* > —1, and 
a*, i € D, satisfy (3.1)-(3.3). By Theorem 3.4 under measure Q the hazard rate 
functions 7 ® are defined by (3.8). 

The family of the equivalent martingale measures for X can be disclosed pre¬ 
cisely. 

Theorem 4.1. Measure Q is the martingale measure for process X if and only if 

c i (t)+a i (t)a*(t)+'yf(t)hi(t)=0, t>0,ieD. (4.1) 

Proof. This result is well known, see e.g. Bellami and Jeanblanc (2000), Proposition 
3.1. The proof follows from Theorem 3.1 and Theorem 3.4. 

Let measure Q be defined by (3.5)-(3.6). Then, by Theorem 3.4, part (c), under 
measure Q the process 

X(t) = T c+aa * (t) + N h {t) + W a (f), t > 0, 

is again the jump-diffusion process with switching regime. The martingale condition 
(2.23) of Theorem 3.1 becomes (4.1). 

The theorem is proved. □ 

The relative entropy functions Hi(t ), t > 0, i £ d, of the martingale measure 
Q w.r.t. IP solve system (3.13) with functions at specified by (3.14) and (3.15). 
Driving parameters c*,h*,a* and switching intensities yf,7p, i G D, satisfy (3.8) 
and (4.1). 

Consider the following examples when the equivalent martingale measure Q is 
unique. 

Example 4.2 (Jump-telegraph process). Consider process X missing the diffu¬ 
sion component, 

X(t) = T c (t) + N h {t), t> 0, 

see (2.6)-(2.7). 
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Assume, that hi{t) ^ 0 almost everywhere, and hi(t) is of the opposite sign 
with d(t) : 


d(t)/hi(t) <0, t> 0, 

i £ D. 

(4.2) 

Moreover, let functions d/hi, i £ D, be locally 

integrable and 


cAu ) 


(4.3) 

/ 1 (\ du ~ °°’ 

J 0 h i\U) 

i £ D. 


Then, by Theorem 4.1 the equivalent martingale measure Q exists and it is unique 
with the hazard rate functions of interarrival times defined by 

yf (i) = - Ci(t)/hi(t ) >0, t > 0, i £ D. (4.4) 


Here (4.3) is the non-exploding condition for measure Q. The corresponding mea¬ 
sure transformation is determined by the functions c*, h*, i £ D, 

C* ( t ) = yf (t) - 7 ?(t), h*(t) = -1 + yf (i)/yf(*), t > 0, i £ D, 

if yf(t) > 0 a.e., see (3.8) and (4.4). * 2 

The entropy functions Hi(t), i £ D, solve system (3.13) with 


ai(t) 



7?0) + 7?(«) In 


7 ? 7 ) 

7 ?(«) 


(u)du, 


where yf are defined by (4.4). The survival functions Ff, i £ D, are defined in 

(2.13). 

If the inequalities (4.2) do not hold, the martingale measure does not exist. 

In particular, if the parameters hi ^ 0, are constant such that d/hi < 

0 , i £ D, with exponentially distributed inter-switching times, yf = Aj, i £ D : 
then by (4.4) under the martingale measure Q the inter-switching times are again 
exponentially distributed with the switching intensities A'® = — d/hi , i £ D. If 
d = 2, then the closed form of the entropy functions is found, see Remark 3.6, 
formulae (3.17)-(3.18), and more detailed analysis in Section 5 below. In this case 
the unique martingale measure Q is defined by the Radon-Nikodym density (3.5)- 
(3.6) with constant c* and h*: 

c* = Xi~X *, h* = -1 + ^, ( 4 . 5 ) 

'n 

where A* = — d/hi >0, i £ D = {1,2}, are the new alternating intensities of the 
inter-switching times, see e.g. Kolesnik and Ratanov (2013). 

Formulae (3.17)-(3.18) for the entropy functions Hi{t), t > 0, hold with 


b,.= 




i £ D = {1,2}. 


4f, in contrary, hi(u) = 0 on a whole interval, u £ (a. h) C [0, oo), then, due to (4.1) with 
= 0 , we have no martingale measures (if Ci ^ 0 on the interval (a, b)), or infinitely many ones 
(if Ci = 0 with free fragment of hazard rate function 7 ^). 

2 Measures Q and P are equivalent. If P-distribution of the interarrival times has a “dead” 
zone: 7 f(u) = 0, u E {a,b), for some time interval ( a,b) C [0,oo), then due to (3.8) and (4.1) 
for any martingale measure Q the hazard rate function 7 ^ also vanishes on (a, 6 ), 7 ®(u) = 0 and 
C{(u) = c*(u) = 0, u G ( a,b ). 
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Example 4.3 (Diffusion process). Consider the diffusion process missing the jump 
component and switching, c is locally integrable and functions <r, c/cr are locally 
square integrable. Assume that cr(u) 0 a.e. 

Let Q be an equivalent measure. In this case the Radon-Nikodym derivative of 
Q is defined by 

= exp ^ a*(u)dB{u) ~ \ cr*(u) 2 du^ , 


where cr* is the locally square integrable function. By Girsanov’s Theorem the 
process 

B = B— f a* (u)du 
Jo 

is Q-Brownian motion. Hence, process X takes the form 

X(t) = I [c(u) + cr(u)a*(u)]du + f a*(u)dB(u), t> 0. 

Jo Jo 

This is a martingale if and only if a a* = — c. 

By (3.11) the relative entropy H(t ) of Q w.r.t. P is 


H{t) = E c 


In -—(f) 
dP v ’ 


=Er 


cr*(u)dB(u) — — I cr*(u) 2 du 


Uo 


=Er 


a* (u)dB(u) + — cr*(u) 2 du 


L^o 

t 


f <r*(u) 2 du. 

Jo 


Therefore, the relative entropy of the (unique) martingale measure is given by 



t > 0. 


Remark 4.4 (Diffusion process with switching tendencies and diffusion 
coefficients). Consider the case of the diffusion process, 


X(t) = T c (t) + W a (t), t> 0, 

with the switching tendencies q = Cj(f) and diffusion coefficients <7* = cr,(t) ^ 
0 , t > 0 , * £ D , where the jump component is missing. 

In this case there are infinitely many equivalent martingale measures. 

Theorem 4.1 shows that the measure transformation defined by \j r t = £t{X*), 

see (3.5) with h* = 0, c* = 0 and a*(t) = —d{t)/<Ji(t ), t > 0, i £ D, eliminates 
the drift component similarly as in Example 4.3. Under measure Q process X 
becomes the martingale of the form X(t) = W a (t), t > 0, whereas by (4.1) the 
inter-switching times are arbitrary distributed. Here W a is the stochastic integral 
(2.8) based on Q-Brownian motion B. 

This model has been analysed in Elliott et al (2005) by using the Esscher trans¬ 
form under switching regimes. This transformation does not affect the distribution 
of inter-switching times and the corresponding equivalent martingale measure is of 
the minimal relative entropy, see Elliott et al (2005), Proposition 3.1. 
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In the next section we study in detail the jump-diffusion model with switching 
regimes based on the Markov underlying process e. We have discovered that in this 
case the Esscher transform does not produce the minimal relative entropy. 


5. Esscher transform and minimal entropy martingale measure 


Typically, the jump-diffusion model with switching regime has no martingale 
measure or it has infinitely many. The rare examples of the unique martingale 
measure are presented above (Example 4.2 and Example 4.3). In this section we 
discuss the case when the infinitely many martingale measures exist and discuss 
some methods to select one. The first method is based on the so-called the Esscher 
transform under switching regimes. 

Let X = T c (t)+N h (t)+W a (t), t > 0, be the jump-diffusion process with switch¬ 
ing regime, see (2.6)-(2.8), defined on the filtered probability space (f2, T, F t , IP)- 

Let dj / 0, i G D, a. s. The case with missed diffusion (ct, = 0, i G D) is 
analysed in Example 4.2. 

To choose an equivalent martingale measure by a reasonable way consider the 
deterministic measurable functions di = 8i(t), t > 0, i G D, which define the 
regime switching processes dj = dj(t), t > 0, i G D, similarly as in (2.1). Let 
measure Qe (equivalent to P) be defined by the density 


exp 

(j> (s)d y( s ) 


Ep 

exp ( 

J> (s )dF( s )) 

i n 



(5.1) 


Here Y(t) = T c ~ a ^ 2 {t ) + N l < 1+h \t) + W a (t), t > 0, see (2.10), and JF e t is the 
P-augmentation of the natural filtration generated by e. This particular choice of 
the new measure is named a regime switching Esscher transform (or exponential 
tilting), see Elliott et al Elliott et al (2005). 

It is easy to see that 


Ero 


exp ( / #t(s)dY(s) 




exp (T e ^ c - a2 / 2 \t) + N 8 ln ( 1+ft -) (t)j exp ('t 


e 2 a 2 /2 


Therefore the Radon-Nikodym derivative of the Esscher transforms is given by 



\^ t = exp ( W 8 °(t)-T e 2 ° 2 / 2 (t )) 


which corresponds to Radon-Nikodym derivative (3.5) with 


= dr 


= 0, h* = 0, 


i G D. 


(5.2) 


Observe, that by Girsanov’s Theorem (see Theorem 3.4, equation (3.8)) due to 
(5.2) the distribution of inter-switching times are not changed under such defined 
measure, = 7 P . 

Hence, due to (5.2), the martingale condition (see Theorem 4.1, equation (4.1)) 
can be written as 


0i(t) = 


Ci(t) + 'yf(t)h i (t) 

<Ji(t) 2 


t> 0, i G D. 


(5.3) 


It is known that the Esscher measure transform defined by (5.1) with parameters 
di, i G D, determined by (5.3) corresponds to the minimal relative entropy , see 
Elliott et al (2005), Proposition 3.1. The similar approach with the Esscher measure 











ON JUMP-DIFFUSION PROCESSES WITH REGIME SWITCHING 


17 


transform produces the minimal relative entropy in the case of Levy processes (see 
Esche and Schweizer (2005); Fujiwara and Miyahara (2003)). 

For our model based on Brownian motion with jumps and with switching regimes 
the Esscher transform does not produce the minimal relative entropy. 

In the rest of this section for the sake of simplicity, we consider the Markov 
case with d = 2 , when the alternating distributions of inter-switching times are 
exponential both under measure P and under an equivalent measure Q, i.e. 

qf = Aj = const >0, 7 ? = A* = const >0, i £ D = { 1, 2}, 

and the driving parameters , i £ D = {1,2}, are constant. Here measure 

Q is defined by (3.4)-(3.6) with constant parameters c*,h*,a*, i £ D = {1,2}, 
satisfying (3.1)-(3. 3). 

To analyse the set of equivalent martingale measures from the viewpoint of the 
relative entropy we are looking for the solution of the integral equations (3.13). 
Since this jump-diffusion process X is bounded, by Theorem 2.1 of Fritelli (2000) 
there exists a unique minimal entropy martingale measure. 

As it is shown in Remark 3.6, in this case the relative entropy functions are 
defined by (3.18): 


H 1 (t)=H 1 {t-,X* 1 ,X* 2 ) =Bt + A 1 
H 2 (t)=H 2 (f,Xt,X* 2 ) =Bt + A 2 


\ _ g — (Aj+Aj)t 
\ _ q — (M+^Sl ) 4 


t > 0 , 


where A±, A 2 and B are defined by (3.17) and (3.15). 


Remark 5.1. Note that B = 0 (or, equivalently, b± = b 2 = 0) if and only if process X 
is already a P-martingale. Indeed, b\ = b 2 = 0 if and only if a* = 0, i £ D = {1, 2}, 
and A* = A*, i £ D = {1,2}, (see (3.15)) , or equivalently, Q = P. Hence 
A 1 = A 2 = 0 and H^t) = H 2 (t) = 0. 


Remark 5.2. Let the jump-diffusion process X be a Levy process, i.e. the alterna¬ 
tion is missing and ci = C 2 = c, hi = h 2 = h 7 ! 0 , o\ = a 2 = a 7 ^ 0 are constant. 
This is the case of a Markov jump-diffusion process. 

Let the new measure Q be defined by (3.4)-(3.6). 

Therefore, by (3.18) the relative entropy functions are identical and linear in t, 

Hi ~ H 2 — Bt 1 t ^ 0, 

where, due to (3.17), B = b = X — X* + X* lnA*/A + (ct*) 2 /2, and Ai = A 2 = 0. 
Here A = j F and A* = 7® are the constant jump intensities under measure P and 
measure Q respectively; the parameter <7* satisfies martingale condition (4.1): 

* c + A *h 

a = -. 

(7 

In this case the martingale measure with the minimal relative entropy is defined 
by the jump intensity A* which satisfies the algebraic equation: 

&'(A*) = InA*/A + —2 (t + A*) = 0. (5.4) 

The latter equation is equivalent to 

c + (3* a 2 + Xhexp(/3*h) = 0, 


(5.5) 
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where the following change of variables A* = A exp (j3*h) is applied. In this par¬ 
ticular case of Levy process X, equation (5.5) coincides with condition (C) of Fu- 
jiwara and Miyahara (2003), which gives the minimal relative entropy H(t) under 
a measure defined by the Esscher transformation. In this example equation (3.18) 
becomes 


Hit) 


r 2 

A(1 - exp(/3*/i) + (5*h exp(/J* h)) + ^ + Aexp(/3 *h)^ 


(5.6) 


where /?* is the (unique) solution of (5.5). Equation (5.6) corresponds to equation 
(3.9) from Fujiwara and Miyahara (2003). 

In the case of the jump-diffusion process with alternating parameters, such coin¬ 
cidence is not available. In this case the minimal entropy functions depend on the 
initial state and they have a bit more complicated behaviour. 


Observe that functions b\ = 6i(A^, cr*), b 2 = b 2 ( A 2 ,u 2 ) are expressed by sum¬ 
ming up of the two nonnegative and convex functions, /( x) = ci — x+x ln(;c/a), x > 
0, (with A* for x) and g{y) = y 2 /2 (with a* for y ). Hence, bi and b 2 are nonneg¬ 


ative and convex. Therefore, function B = 1?(A}, A 2 , <j*, er 2 ) = 


A 2 b\ + A \b 2 
A} + A 2 


is also 


nonnegative. 

We analyse the relative entropy functions H\ = Hi{t) and H 2 = H 2 {t) for 
small and big times t separately. These functions possesses the following time- 
asymptotics. 


Proposition 5.3. Let the relative entropy functions H i and H 2 be defined by 
(3.18). Thus, 

Hi{t) ~ bit, H 2 (t) ~ b 2 t, t -)> 0, (5.7) 


and 

Hift) ~ Bt + Ai, iJ 2 (t) ~ Bt + A 2 , 
Proof. As t —> 0 formulae (3.17)-(3.18) lead to (5.7): 


oo. 


Hi (t) — Bt + Ai 


1 — e 


— (AJ-l-AJ)t 


[B + A^(AJ + A 2 )] t = bp, 


(5.8) 
i&D = {1,2}. 


Long-term asymptotic (5.8) is evident. 


□ 


We choose the equivalent measure minimising the relative entropy function H(f) 
under the martingale condition (4.1), i.e. A*, a*, i £ D = {1,2}, satisfy the 
following relations: 

j c\ + A* • hi + c* • <ti = 0, 
c 2 + A 2 • h- 2 + u 2 • <t 2 = 0. ^ ^ 

If measure Q is constructed by way of minimising b\ and b 2 , we say that Q is 
the short-term minimal entropy martingale measure (MEMM), see (5.7); measure 
Q is the long-term MEMM, if Q minimises B, see (5.8). 

The short-term and the long-term minimal entropy equivalent martingale mea¬ 
sures are defined as the solutions of the following optimisation problems subject to 
martingale condition (5.9): 

• the short-term MEMM is defined by solving the problem w.r.t. a* and 

A*, i & D = {1,2}, 


6i = Ai — A} + A} In A{/Ai + K) 2 /2 -> min, 
b 2 = A 2 — A 2 + A 2 lnA 2 /A 2 + (<t 2 ) 2 /2 —> min. 


(5.10) 
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the long-term MEMM is defined by solving the problem w.r.t. a* and 
X*, i €£> = {1,2}, 


B = 


A^l ^1^2 
At+A* 


(5.11) 


Theorem 5.4. The solutions of problems (5.10) and (5.11) subject to condition 
(5.9) exist and they are unique. 

Proof. If Ui 7 ^ 0, then (5.9) gives 


Ci + A* hi 


(5.12) 


Hence problem (5.10) is equivalent to minimisation of the function 

b t = bi( A*) := Xi - X* + X* ln(X*/Xi) + ±- 2 (a + X*h,) 2 , i G D = {1,2}. 

If, additionally, hi = 0, then problem (5.10) has the unique solution A* = Aj. In 
this case we return to the Markov modulated diffusion process, see Remark 4.4 and 
Elliott et al (2005). This confirms again that in this case the minimum of entropy 
and the Esscher transform (5.1) lead to the same result. 

If <Ji = 0, hi ^ 0, then we have jump-telegraph model (see Example 4.2). In 
this case the martingale condition (5.9) gives A* = — Ci/hi , which corresponds to 
the unique equivalent martingale measure. 

On the contrary, if < 7 * 7 ^ 0, hi 7 ^ 0, then the minimal entropy and the Esscher 
transform (5.1) give different results. Observe that bi, i G D = {1,2}, with a* 
given by (5.12) can be rewritten as 

b t = bi (A*) := Xi - X* + X* In (A*/A,) + (A* - a if , * G D = {1,2}, (5.13) 

where Cf = hf/af > 0 and a* = — Ci/hi >0, isD = {l, 2 }. 

Differentiating (5.13) we have: 


&i(A{) =— \(c.\!h\ + A}) + ln(A{/Ai), 

a i 

^ 2 (^ 2 ) = — f (C 2//12 + X 2 ) + ln(A 2 /A 2 ). 


(5.14) 


We remark that functions b[ and b' 2 vary monotonically from —00 to +00 as A} 
and X 2 increase from 0 to + 00 . Hence, the system b\ (A{) = 0, b 2 (X 2 ) = 0 (and 
minimisation problem (5.10)) has the unique solution. 

Moreover, if A* < — a/hi , then the solution A* of the corresponding equation 
&'( A*) = 0 satisfies A*, X* > X t , i G D = {1,2}; if A* > —Ci/hi, then A* < A*, i G 
D = {1,2}. Therefore, the solution A* of (5.10) is always between Aj and — Ci/hi , 
i G D = {1,2}, whereas the Esscher transform gives A* = A„ i S ff = {1, 2}. 

We solve the problem (5.11) by differentiating again. System 


dB 

SAf 


= 0 , 


dB 

<9AS 


= 0 


is equivalent to 


$1 := (AJ + Al)&i(AD + 6 2 (A*) - h(Xl) =0, 
$ 2 := (AI + X* 2 )b' 2 (X* 2 ) + AJ) - b 2 iX* 2 ) =0. 


(5.15) 

(5.16) 
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By a similar reasoning as before, one can easily see that system (5.15)-(5.16) also 
has (unique) solution. Indeed, owing to (5.13) and (5.14) we have 

= A; ln^/Ar) + b 2 { AS) - X 1 T A} + ^(AJ - m) 2 + C 2 (A{ - ai )(AS T m). 
Thus, 

^ = 1 + ^ + C 1 2 (At + A^) > °. 

Hence, function $i increases monotonically from —oo to +oo as A} increases from 
0 to Too, which means that equation (5.15) has the unique solution A} = <j)(X 2 ) > 0 
for any fixed positive AS- Therefore system (5.15)-(5.16) is equivalent to 

b[(cf>(\* 2 )) + b' 2 (\* 2 )=0. 

Differentiating this identity and (5.14) one can see that 

C'2 T 1/AS 


4 >'(AS) = -- 


< 0. 


Cl T 1/0(AJ) 

Hence, function &i(0(^2)) decreases as AS goes from 0 to Too, whereas b' 2 (X 2 ) 
strictly increases from —oo to Too. Therefore system (5.15)-(5.16) has the unique 
solution. 

Then, the second derivatives of B are 


B n =- 
B 22 = 
B\2 = B 21 = 


AS 


AT T A. 


-b'{( AI)- 


2AS 


XI)- 


A^ T AS 

At 


(At T A^) 3 
2At 

(AJ T AS) 3 


r $i(At,AS)T 


$i(At,AS), 

$ 2 (At,AS), 

$ 2 (At,AS). 


A: 


d 2 B 

Here B - = Wdx 


(AtTAS ) 3 v 11 2/ (AtTAS ) 3 

*, i,j e D = {1,2}. Hence, point (AS, AS), fitting for system 


(5.15)-(5.16), gives the minimum of B. 


□ 


Remark 5.5. Theorem 5.4 and Proposition 5.3 show that the minimal entropy mar¬ 
tingale measure differs from the measure supplied by Esscher transformation. This 
is confirmed by the plots presented in Fig. 5.1 and Fig. 5.2, where an asymmetric sit¬ 
uation is considered. Surprisingly, the minimal entropy is supplied by A}, AS which 
depend on time horizon. Moreover, the entropy functions iJ;(t), t > 0, i £ {1,2}, 
are not linear (cf. (5.6), Remark 5.2). 

In Elliott et al (2007) the Esscher transform is applied to option pricing under 
a Markov-modulated jump-diffusion model. 

In the symmetric case, 

Ai = A 2 , er 1 = cr 2 and Ci = —c 2 , hi = —h 2 , C\h 2 = c 2 hi , (5.17) 

the solution of the minimal entropy problem is constant. 

Proposition 5.6. In the symmetric case (5.17) the problem 

Hi(t; A}, AS) -T min, 

H 2 (t\ A}, AS) min 


(5.18) 
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H,m 



H 2 m 



Figure 5.1. Plot of Hi(t)/t, t > 0, where Hi(t ), i £ {1,2}, is the 
minimal entropy (for Ai = A2 = 1 , ay = <72 = 1, ci = — 1, C2 = 3, 
hi = 1 , h 2 = — 0 . 1 ). 


subject to condition (5.9) has the unique solution A} = A(5 = const, which does not 
depend on time t. 


Proof. By Taylor representation we have 

Hi =bit + Xl{bi - b 2 )t 2 ■ <\> ((AI + A * 2 )t ), 
H 2 =b 2 t + A * 2 (b 2 - bi)t 2 ■ </> ((AI + A * 2 )t ), 


where 


<t>{x) = 

n —0 


(—l) ra+ 1 a;" 
(n + 2 )! 


Notice that in the symmetric case (5.17) function bi and b 2 (see (5.13)) are equal: 
61 (A) = 62 (A) = 6 , A > 0. Hence, the problem (5.18) has the unique constant 
solution, A} = AJ = A*, which corresponds to problem (5.10). The minimal entropy 
is linear, Hi(t ) = H 2 (t ) = 6 t, where 6 = 6 (A*). □ 
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